6 


daa 
NACA TN 3317 LS 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


TECHNICAL NOTE 3317 


DESIGN CONSIDERATIONS FOR WINGS HAVING 
MINIMUM DRAG DUE TO LIFT 
By Warren А. Tucker 


Langley Aeronautical Laboratory 
Langley Field, Va. 


Washington 
December 1954 


ийїї 


v 


TECH LIBRARY KAFB, NM 


ШИ 


NATIONAL ADVISORY СОММТТТЕЕ FOR AERONAUTICS 029 


TECHNICAL NOTE 5517 


DESIGN CONSIDERATIONS FOR WINGS HAVING 
MINIMUM DRAG DUE TO LIFT 


By Warren A. Tucker 
SUMMARY 


The problem of increasing the range of supersonic aircraft by the 
use of twisted and cambered wings is considered, primarily for the pur- 
pose of developing a rational method for the selection of а, design lift 
coefficient. Relations and curves are presented from which a suitable 
selection may be made, depending on the relative importance of maximum 
range and top speed. 


INTRODUCTION 


The inherently low values of lift-drag ratio and the resultant 
short ranges characteristic of conventional aircraft configurations 
operated at supersonic speeds have stimulated research on more efficient 
shapes for use at these speeds. This research has of necessity been 
directed along two main lines, one the reduction of the drag at zero 
lift and the other the reduction of the drag due to lift. Along the 
latter line, one of the most promising developments has been the use 
of twist and camber to reduce the wing drag due to lift for a given 
plan form. Coptributions to this field have been made by Robert T. 
Jones (refs. 1 to 3), E. W. Graham and his coworkers (refs. 4 and 5), 
and S. H. Tsien (ref. 6), among others. 


One of the principal problems with which the reports just mentioned 
are concerned may be stated thus: Given а wing plan form operating at a 
given lift and a given Mach number, find the shape (as expressed by the 
angle of attack, the twist, and the camber) which will result in the 
lowest drag. Attention is directed to the words "at a given lift and a 
given Mach number since the optimum shape will change both with the 
lift and with the Mach number. For апу particular aircraft, however, 
practical considerations will usually dictate the use of a wing having 
fixed twist and camber, even though the aircraft may be required to 
operate over & range of 142+ coefficient as well as a range of Mach 
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number. Out of this situation arises the problem of how best to select 
the one particular shape that will represent the best compromise for all 
flight conditions, which amounts to selecting the lift coefficient ал 
which the drag is to be minimized at a given Mach number. The present 
paper offers some suggestions concerning this selection process. 


SYMBOIS 
Ср drag coefficient, ог. 200 
Dynamic pressure Х Reference area 
Ср, ғ drag coefficient of flat wing at zero lift 
2 
Сту drag coefficient of particular minimum-drag wing at 
O,w 
zero lift 
CDa drag coefficient of particular minimum-drag wing at 
design lift coefficient 
Cp, lowest drag coefficient for a particular minimum-drag 
,W wi 2 | 
ng 
ACh drag coefficient measured from the flat-wing zero-lift 
drag coefficient, Cp - Cp , ‘(used with the various 
2 
subscripts ) 
Ст, lift coefficient, э eee. eee eee 
Dynamic pressure X Reference area 
Ста - design lift coefficient D | ын 
“Те lift coefficient at which drag of particular minimum-drag 
wing equals drag of flat wing 
ET lift coefficient at which CDn, w is obtained 
CLopt,f lift coefficient at which (L/D)nax, f is obtained 
СЕ а lift coefficient at which (L/D)max,~ is obtained 


1/2/ 
Стор - lift coefficient &t which 6 — is obtained 
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(L/D nase, £ 
СЕР max, f 
(L/D max jw 
(13/2 Jo урт 


Kr 


СТе, = дст /до 


а, 


Superscripts: 


*3* 


maximum liftedrag ratio for flat wing 


maximum (11#+)1/2-агав ratio for flat wing 


maximum lift-drag ratio for particular minimum-drag wing 


maximum (1ift)l/?-àrag ratio for particular minimm-dr 
wing 


drag-rise factor for flat wing, ӚСр/ӘС12 


drag-rise factor for family of minimum-drag wings, 
3cp/acr? 


Mach number 


cotangent of sweepback angle of leading edge 


angle of attack 


(L/D) ах, у 
(L/D аах, ғ 


(0/5 = 
11/2 h 


condition of maximum 


condition of maximum 


max, Ё 
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ANALYSIS AND DISCUSSION 


Factors Influencing the Range 


It is anticipated that aircraft operating at supersonic speeds will 
usually be powered by a reaction propulsion system (turbojet, ram jet, 
or rocket, for example), for which the fuel rate (weight of fuel used 
per unit time) will be closely proportional to the thrust, rather than 
to the horsepower, as in the case of-a propeller-driven airplane. For 
such aircraft, the maximum range will be attained by flying at the con- 
dition of (L/D)me, if the flight Mach number is specified, or at the 


condition of (ef) if the ratio of the atmospheric pressure to 


the wing loading is specified (specifying this ratio is roughly equivalent 
to specifying the altitude). These relations are illustrated in figure 1. 


Tbe characteristics of the propulsion system enter into the pre- 
ceding statements only to the extent that the fuel rate is assumed to 
be proportional to the thrust (for example, the changes in available 
thrust with changes in altitude or Mach number do not affect the con- 
siderations). In the case of any specific aircraft, however, the engine 
characteristics will be a primary factor in the determinstion of the 
actual operating point for maximum range. Although a detailed discus- 
sion of this problem is not the purpose of this paper, it is mentioned 
so ав to emphasize the fact that the operating point for a particular 


aircraft may be neither at (L/D)max nor at 11/2 fp max? Put at some 


point between these two conditions. Therefore, after the main assump- 
tions and relations have been stated, the following analysis is divided 
into two main parts, one dealing with (L/D)gg, and the other dealing 


with (1/2/) so that the limits of probable design conditions are 


covered. 


Basic Assumptions and Relations 
The drag coefficient of the flat wing is assumed to be expressible 


as a parabolic function of the lift coefficient for a given Mach number 
and wing plan form: 


Ор = Cro + KeCy* | (1) 
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This equation is plotted in figure 2. As long as the relation between 
Cp and Ст is parabolic, the curve can be regarded as applying either 
to a complete aircraft configuration or to an isolated wing. For con- 
venience, however, the term flat wing is used herein when referring 

to the curve. The drag of the family of wings which at any value of Ст, 
has the lowest possible drag is assumed to be known (from the work of 
ref. 6, for example, or from some other suitable source) and to be given 
by the following equation: 


= 2 
Ср = Ср e + Kyl (2) 
This equation is also plotted in figure 2. 


The first step is to write the equation for the drag of one of the 
family of minimum-drag wings, that is, of a wing which has a fixed twist 
and camber. Мом the addition of small amounts of twist and camber to a 
flat wing introduces increments of drag and lift which do not vary with 
angle of attack. The drag of a wing with fixed twist and camber can 
therefore be represented by the following equation: 


Ср = Cp, y + Ке(СІ - a | (5) 


Where С and С are shown in figure 2. If equation (3) is to 
Од, Im,w 


express the drag of one of the family of minimum-drag wings, the value 
of Cp given by equation (3) must first be equated to the value given | 


by equation (2) for some particular lift coefficient Ста; which is 
called the design lift coefficient. In addition, since the drag given 
by equation (2) is the lowest possible at any value of Ст, then the 


curves represented by equations (2) and (3) must be made tangent at 
Cy, = Ста“ These two conditions serve to determine CDn, w ала Clan yw 


in equation (3), and the equation for the drag of the particular 
minimum-drag wing becomes 


2 
Cp = CDo ғ t ( us Эрэлч + Ке|Ст,- ( - JaA (в) 


Тре relations just discussed are illustrated in figure 2. The discus- 
sion of the curve represented by equation (4), in particular its ге1а- 
tion to equation (1), forms the remainder of the paper. Some miscel- 
laneous relations for the various coefficients are presented in 
appendix A. 
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Design Considerations Based on (1/9), 


Maximum lift-drag ratio.- As mentioned previously, twist and camber 
may be used to produce greater values of the lift-drag ratio, and in 
particular of the maximum lift-drag ratio (Шр) as w than сап be 


3 
obtained from the flat wing. From equations (4) and (1) the ratio of the 
maximum lift-drag ratio for any particular minimum-drag wing (L/D) max yw 


to the corresponding value for the flat wing (L/D) max, ғ can be obtained 
as a function of the two quantities K,/Kp and Ста | Lot р› where 

: | 
CLopt,f 18 the lift coefficient corresponding to (L/D max, р. The ratio 
is as follows: 


К w a O O (5) 
(L/D)max , f 


оценет не 
2|Кебр + к 
Сро, ғ 2 
Vou - Ke 


The ratio given by equation (5) is plotted in figure 5(а). Of some 
practical interest is the fact that for a given value of К/К (which 


amounts to a given plan form and Mach number) there is a particular design 
lift coefficient at which the greatest increase in (L/D) max is obtained. 


(L/D) ах W 

ТГ. ап asterisk is used to denote the condition of maximum ; 
(1/0) ах f 

then the following relation сап be written: 

X 
(L/D) С 
L 
max ,W яа? 
“=. (7) 


(L/P дек, же үк, /Ке 
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Equation (T) is shown as the dashed line in figure 3(a). One criterion 
is thus provided for the selection of Cigs if the only item of concern 


in the particular design under consideration is the maximum possible 
increase in (L/D) to the exclusion of other items (such as increases 


in minimum drag). 


max? 


The numerical values of the (L/D), ratio in figure 5(а) are also 
of interest. As shown in appendix B, values of К,/Кр in the neighbor- 


hood of 0.5 or 0.6 can probably be expected for reasonable plan forms and 
Mach numbers. From figure 3(a), increases in (L/D)ma, of 30 to 40 per- 


cent can therefore probably be realized by the introduction of the proper 
twist and camber. 


Other characteristics.- In connection with the lift-drag ratio, the 
lift coefficient corresponding to (L/D),,, is of importance. The ratio 


of Торо yw to CLopt ,f is given by the following equation: 


2 
CLopt ;W =| eed (8) 
CLopt f Оїорь,г 


This ratio is, of course, always greater than unity. 


In the case of configurations for which the maximum speed capabili- 
ties, as well as the maximum range, are important, the increase in mini- 
mum drag caused by the addition of twist and camber to the flat wing must 
be considered. (Although a detailed study should include the effect of 
Mach number, the change in minimum drag with Mach number is small enough 
so that for the purposes of this paper the constant Mach number case can 
be considered.) This increase in minimum drag can be expressed as 


follows: 
2 
“Dm yw = zb 2 ee) (9) 


c 
Lopt,f 


and the lift coefficient at which the minimum drag occurs can be deter- 
mined from 


С | C 
и. е _ = (10) 
CLopt,f ne CLopt,f 
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If, as before, an asterisk is used to denote the condition of 


L 
maximum Р) пах, „ then equations (8) to (10) become 
{L/D) ах f 
(ас) 2 “Лаг шоо (118) 
CLopt ,f Горе к/к» 

Cp, i 
- 22122 - +5 Bs z (11b) 
Do,f (Caf 1g е) 


(The reason for the seemingly strange result given by this equation 
for the impractical (or ideal) case of Ку —3 О can be understood if a 


plot is made similar to fig. 3(c) but with CL, [Сто as the abscissa. 
р,” 
The case of Kp —>® can be visualized easily by making a sketch similar 
to fig. 2.) | 
: Ко 
С | 
va 1 = = (lic) 


Equations (8) to (11) are plotted in figures 3(b) to 5(4). It is instruc- 
tive to compare the maximum increase in (L/D), with the corresponding 
increase in minimum drag, that is, the values given by the dashed lines 
of figures 5(а) and 3(c), which come from equations (7) and (11b). These 
values arè plotted in figure 4 and show that for a realistic range of 

Ky / Ke (greater than 0.4) the percentage increase in minimum drag is 


always greater than that in (L/D) лу if the design condition is chosen | 
to give the greatest possible increase in (L/D),4,,. This situation is 


very undesirable, of course, if the top speed of the configuration under 
study is important, since the top speed is greatly dependent on the mini- 
mum drag. 


Evidently if the top speed of the configuration is important, a more 
satisfactory design condition would be one which would yield something 
less than the maximum possible (L/D), increase but which would keep 
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the minimum-drag increase within acceptable limits. Such a design соп- 
dition is possible, as shown by the data presented in figure 5. As an 
example, consider a curve with K,/Kp = 0.5. The maximum possible 


(L/D) ах increase of 41 percent is achieved only at the expense of a 


50-percent increase in minimum drag. However, if the design condition 
is modified so that а (say) 30-percent increase in (L/D)y,, is obtained 


(by choosing a value of Cz, | їсрь,е of about 0.68), then the minimum- 


drag penalty is decreased to 12 percent. Similar considerations apply, 
of course, to any of the other curves with К,/Кр constant. Ав indi- 


cated in the example just discussed, the lines of constant Ora /Сторь Р 
2 


give the proper value of design lift coefficient to use once a satis- 
factory design condition, as measured by the gain in (L/D) ах balanced 


against the penalty in minimm drag, is selected. 


1 
Design Considerations Based оп (2/2) ах 


Considerations analogous to those developed in the preceding section 


entitled "Design Considerations Based on (L D) ах Can be expressed on 


the basis of (11/2 jp). ах: This development proceeds in a manner similar 
to that for the case of ( L/D) max: 


Maximum value of 11/2) ‚- The equation for (2/2) согге- 


sponding to equation (5) is 


C M CE 


(22) ur (C fon) ase 


| 1 (12) 
ARS TET REP BES 
For the flat wing, the following relations hold: 
a = ПЕ: И "EE. Ж 
так ae ep, De рые (15) 


c _ | 20,2 
Lopt,f Ke 
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It is emphasized that. CLopt " has the same meaning here that it had in 
2 
the previous development; that is, it is the lift coefficient corresponding 


to (L/D) аху ( not (о) ах +): 


The ratio given by equation (12) is plotted in figure 6(a). Ifa 

(11 /2/ Dess 
J 

(17/7 /p) 


max ,f 


Gouble asterisk is used to denote the condition of maximum 
then the following relation holds: 


(11/2 f) ах | 
и opet (к) 


Equation (14) is plotted as the dashed line in figure 6(a) and provides 
a criterion for the selection of Ста? the basis of selection in this 


case being the development of the maximum possible increase in 
(т!/ а Го) ы to the exclusion of other factors. 


(14) 


Other characteristics.- The lift coefficient corresponding to 


1/2 1 1 
(2 №) ах м is designated herein by Сто ж The ratio of Сте уу 
to “1 - is given by the following equat ión, which corresponds to 


ро, 
equation (8) for the case of (1/2): 


The equations for minimum drag and the lift coefficient at which minimum 
drag occurs (eqs. (9) and (10)) apply to the present case as well as to 
the previous one. 


If, as in equation (14), a double asterisk is used to indicate the 


(11/ 2/D max у 
condition of maximum „еа, then equations (15), (9), апа (10) 


(12/2 /p) 


take the following forms: 


max yf ere : 2. sees 
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opt ,w CL _ 1 M d: opt ,w (16a) 
L 


С 
Торы г opt ,f үк, /Ke {5 ы; 


= 
дам) 0h ИНН ц : = (16b) 
D 2. 9 Ke 
0,? 9(Ста/Сгорь,г ) 
8 хэс Е 1 - Be 
low = Та ЕНЕ! ЗЕН EY Ке, 4 кы (16c) 


C C CL./C Танга Y. 
Lopt Pa Lopt f ? 14] Lopt yt 15 Ст Гор, Ї 


^ 


Equations (9), (10), (15), and (16) are plotted in figures 6(b) to 6(d). 


1/2 
A comparison of the maximum increase in (т; 1/ №). with the corre- 
sponding increase in CD is given in figure T. The values in figure 7 


come from equations (14) and (16b). The comparison is rather more favor- 
able than the corresponding comparison in figure 4, but in some cases the 


1/2 
balance between (z / №) 
шах, 


be inacceptable. In such cases a more suitable design compromise can be 
selected from the curves of figure 8, which corresponds to figure 5 for 
the (L/D) max case. Considerations similar to those discussed in con- 


nection with figure 5 apply to figure 8. 


and CDn Эн represented by figure Т may 
2 


Comparison of Designs Based оп (L/D),,. and on (т; е) ы 


Because of the existence of at least two design bases (that is, 


(L/D) ee and (11/2 /p) „.„) the question naturally arises as to what 


extent a design which emphasizes one of these factors is penalized with 

respect to the other. Some indication may be had from the answer to the 

following specific question and its converse: If C /Lopt ,f is chosen 
L/D LL [2 D xx 

to give (2/0) аху ‚ what percentage of Ce Pas will be 
(L/D) max ,£ (11/2 /D) ах, ғ 
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realized? The answers to these questions are shown in figure 9, together 
with the minimum-drag increase and the design lift coefficient for each 
design condition. The first point to be made from the figure is that,- 


1/2 ) 
regardless of which factor (that is, (L/D),,, or (1 /D/y,) is 


maximized, the increase in the other is always at least 92 percent of the 
maximum possible increase. The second point is that the minimum-drag 


increase resulting from maximizing (11/2 /p) ах is much less than that 
resulting from maximizing (L/D) max Therefore, a design lift coeffi- 
cient chosen to give the maximum possible increase in (11/2) „„ will 


result in a wing which will develop at least 93 percent of the maximum 


possible increase in (L/D) пах but which has а much smaller minimum drag 


than if the design had been chosen to give the maximum possible increase 
in 65 дэр eee Such a design lift coefficient therefore suggests itself 


ав а good choice as long as the minimum drag is not too critical. For . 
cases in which closer attention must be paid to the minimum drag, a more 
suitable compromise can be selected from figures 5 and 8. 


CONCLUDING REMARKS 


Under the assumption that, for a given wing plan form, the lowest 
drag obtainable by the use of twist and camber is known, information has 
been presented which will allow the rational selection of a design lift 
coefficient. It is suggested that, for a given Mach number and wing plan 
form, a design lift coefficient which gives the largest possible increase 


in the maximum (11¢t)* er ratio is a reasonable choice as long as the 
minimum drag is not too critical. If the increase in minimum drag due to 
twist and camber must be kept very small, а somewhat lower design lift 
coefficient should be used. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Даана ess 
Langley Field, Va., September 8, 1954. 
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APPENDIX A 
MISCELLANEOUS RELATIONS FOR THE VARIOUS COEFFICIENTS 


Certain relations, although not directly bearing on the main purpose 
of the paper, are nevertheless of incidental interest and are easily 
derived. For example, if the values of Cp given by equations (1) 


and (4) are equated and the resultant lift coefficient is called Cr. 
(see fig. 2), the following simple result is obt&ined: 


Сте 1 
Gr, 2 | (Al) 


Other similar relations can be obtained. For example, if ЛС is 


used to denote the drag coefficient measured from the flat-wing zero-lift 
drag coefficient as a reference, that is, 


Абр = Ор - Оро p (A2) 


then these relations are as follows (see fig. 2 for the physical meaning 
of the various coefficients): 


Тим. Anw _ 122 Ку 
Ста АСрд Kp 
AC 
Зи Ру (45) 
АСро у Ке 
= Do Е ( 5 3 Ку 
АСра Kp /| Ќе 


The preceding relations are shown in figure 10. 
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APPENDIX B 
PROBABLE RANGE OF К/К 


In reference 2, Robert T. Jones gave an equation for the minimum 
drag (at a fixed lift) of a slender wing lying near the center of the 
Mach cone. This equation was also given later by Adams and Sears in 
reference 7, with some discussion of the derivation. For wings with _ 
zero tip chord, the equation may be written in the notation of the pres- 
ent paper as follows: 


Lk, = de 1+ 2(pm)~ (в1) 


n Bm 


The drag for the corresponding flat wing, if the leading edge is assumed 
to develop no thrust, is simply 


1 mm m 
p - Bor. 2: (82) 
so that 
Ky ВОТ 1 + 2(ва)” (85) 


Ke A/m npm __ 


Values of BCT can be obtained from various sources, such as 


reference 8. 


Values of К,/Кр calculated from equation (B5) are plotted in fig- 


ure ll. Also included are some more exact values for the conical trian- 
gular wing (A/m = №), taken from reference 6. The agreement between | 
these values and those calculated from equation (B3) is good. Although 
equation (Bl) is not strictly applicable to a wing with A/m = 6, the 
values in figure 11 are included to show the trend of the variation of 
К„/Ке with A/m. The lower values of K,/Kp obtained for the higher 
values of A/m indicate that the use of proper twist and camber should 
be more beneficial for arrow-type plan forms than for diamond-type plan 
forms. Values of K,/Kp in the neighborhood of 0.5 or 0.6 can probably 
be expected for reasonable plan forms and Mach numbers. 


—ытт . 


As mentioned previously, the values of К,/Ке presented in fig- 


ure ll are based on the assumption that no leading-edge thrust is devel- 
oped in the flat wing. For increasing amounts of leading-edge thrust, 
the values will move closer to unity. 
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Figure 1.- Range relations for a typical set of aerodynamic characteristics. 
Fuel rate proportional to thrust (turbojet, for example). 
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Figure 2.- Basic relations. 
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Figure 5.- Variation of several aerodynamic characteristics of minimum-drag wings with design lift 
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Figure 5.- Variation of (L/D),,, ratio with Cp ratio. 
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Figure 6.- Variation of several aerodynamic characteristics of minimm-drag wings with design lift 
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Figure 8.- Variation of G 2] D) ах ratio with Cp ratio. 
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Figure 9.- Comparison of increases in (L/D) max and 11/ јо) ыл» and corresponding increases in Ср, 
for two design conditions. 


LTEC КЕ VOVN 


NET 
y 


NACA TN 3317 


14 


12 


ч 
Су 


Ratios of соей чел (5 
Qo 


Figure 10.- Relations among characteristics of minimum-drag wings 
(eqs. (А1) and (A3)). 
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Figure 11.- Values of K,/Kp for wings with zero tip chord. No leading- 
edge thrust. E 
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